Control and Tomography of a Three Level Superconducting Artificial Atom 
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A number of superconducting qubits, such as the transmon or the phase qubit, have an energy level 
structure with small anharmonicity. This allows for convenient access of higher excited states with 
similar frequencies. However, special care has to be taken to avoid unwanted higher-level populations 
when using short control pulses. Here we demonstrate the preparation of arbitrary three- level 
superposition states using optimal control techniques in a transmon. Performing dispersive read- 
out we extract the populations of all three levels of the qutrit and study the coherence of its excited 
states. Finally we demonstrate full quantum state tomography of the prepared qutrit states and 
evaluate the fidelities of a set of states, finding on average 96%. 

PACS numbers: 42.50.Ct, 42.50.Pq, 78.20.Bh, 85. 25. Am 



Spin 1/2 or equivalent two- level systems are the most 
common computational primitive for quantum informa- 
tion processing pQ. Using physical systems with higher 
dimensional Hilbert spaces instead of qubits has a num- 
ber of potential advantages. They simplify quantum 
gates [2], can naturally simulate physical systems with 
spin greater than 1/2 [3 , improve security in quantum 
key distribution [4j [5] and show stronger violations of lo- 
cal realism when prepared in entangled states [6j[7]. Mul- 
tilevel systems have been successfully realized in photon 
orbital angular momentum states [8, 9 , energy-time en- 
tangled qutrits [10] and polarization states of multiple 
photons [IT] . Multiple levels were used before for pump- 
probe readout of superconducting phase qubits [T2HT4] , 
were observed in the nonlinear scaling of the Rabi fre- 
quency of DC SQUID's [T5HT8] and were explicitly popu- 
lated and used to emulate the dynamics of single spins [3] . 
In solid state devices, the experimental demonstration 
of full quantum state tomography [10 of the generated 
states, i.e. a full characterization of the qutrit, is cur- 
rently actively pursued by a number of groups. 

In this work, we use a transmon-type superconduct- 
ing artificial atom with charging energy E c /2tt — 298 ± 
1 MHz and maximum Josephson energy EJ ia, */2ir = 38 
GHz p~9j [20] embedded in a coplanar microwave res- 
onator of frequency u r /27r = 6.9421 ± 0.0001 GHz in 
an architecture known as circuit quantum electrodynam- 
ics (QED) [HI [22]. In circuit QED, the third level 
has already been used, for instance, in a measurement 
of the Autler-Townes doublet in a pump-probe experi- 
ment [23] [24] . It has also been crucial in the realization 
of the first quantum algorithms in superconducting cir- 
cuits [25] and is used in a number of recent quantum 
optical investigations, e.g. in Ref. [26 . Also, quantum 
state tomography based on dispersive readout [2TJ [27] of 
a two-qubit system has been demonstrated [28] and used 
for the characterization of entangled states [25[ [29] . In 
our realization of three level quantum state tomography, 
we populate excited states using optimal control tech- 



niques [30] and read out these states using tomography 
with high fidelity. We determine all relevant system pa- 
rameters and compare the data to a quantitative model 
of the measurement response. 

The transmon coupled to a single mode of a resonator 
is well described by the linear dispersive hamiltonian p~9j 
|2T] approximated to second order 
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where the transmon transition frequency uoi = oo\ — ooq 
is largely detuned from the resonator. Here, a^) is the 
annihilation (creation) operator for the photon field and 
Xn = 9n/^n- 9n denotes the coupling strength to the 
transmon transition n o n+1 and A n the detuning of the 
same transition from the cavity frequency. We extracted 
go/2ir = 115 =b 1 MHz from a measurement of the vacuum 
Rabi mode splitting [22] . Coupling constants g n of higher 
levels were explicitly determined in time resolved Rabi 
oscillation experiments, where g n = r\ n g§ « y/n + \g§ 
due to the limited anharmonicity of the transmon [T9] . 
For the uo\2 transition we experimentally determined r]i — 
1.43 ± 0.04. Using flux bias, we detune the qubit by 
Ao = cjoi — r — —1.319 ±0.001 GHz from the resonator. 
The non-resonant interaction with the transmon in state 
|n) leads to a dispersive shift s n = —(Xn — Xn-i) m 
the cavity frequency. Measuring the in-phase quadrature 
amplitude (a^) of microwaves transmitted through the 
resonator [Fig. [TJa)] at a chosen detuning A rm = uo r — 
uo m = 5.1 MHz of the measurement frequency uo m from 
o; r , allows to extract the population of the transmon state 
n. 

To prepare arbitrary superposition states of the low- 
est three levels of the transmon we use optimal control 
techniques in which two subsequent DRAG pulses [30] of 
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FIG. 1: (color online) a) Calculated in-phase transmission 
through the resonator for transmon states n = 0,1,2. The 
dashed curve indicates the bare resonator response. The ver- 
tical blue arrow indicates the detuning A rm of the measure- 
ment tone, b) Pulsed I quadrature measurement responses for 
prepared states |0), |1) and |2). c) Measured dispersive shifts 
s n versus transmon transition frequency Cc?ni- The solid lines 
are calculated within the linear dispersive approximation. 

standard deviation 3 ns and total length 12 ns are ap- 
plied to the qubit at the ujqi and UJ12 transitions. We 
extend the technique described in [30] for the two lowest 
levels of the transmon to three levels using quadrature 
compensation and time-dependent phase ramps [31] to 
suppress population leakage to other states and to ob- 
tain well denned phases. 

For a first characterization of the readout of higher 
levels, the transmon is prepared in one of its three low- 
est basis states \i) (i = 0,1,2). After state preparation, 
a coherent microwave tone is applied to the cavity and 
the state dependent transmission amplitude is measured, 
Fig. [ljb) . The amplitude of the tone was adjusted to 
maintain the average population of the cavity well below 
the critical photon number n cr it = A§/4^q = 25 [2T] . 
The time dependent transmission signals are characteris- 
tic for the prepared qubit states and agree well with the 
expected transmission calculated based on Cavity-Bloch 
equations [27]. We have generalized the formalism pre- 
sented in Ref. 27 to three levels to quantitatively model 
the dispersive measurement. From the fits in Fig. [ljb) , 
we have extracted the state dependent cavity frequency 
shifts 5 ,i, 2 /27r = 10.0, 5.9, 3.4±0.1 MHz, which are found 
to be within 0.1 MHz of the values calculated from inde- 
pendently measured hamiltonian parameters. Also, the 
dispersive frequency shifts s n measured in this way agree 
well with the linear dispersive model over a wide range 
of transmon transition frequencies cjqi, see Fig. [TJc). 

The frequency shifts can also be obtained more di- 
rectly by measuring the transmission amplitude over a 
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FIG. 2: (color online) a) Measured Q quadrature of the res- 
onator transmission versus time and measurement detuning 
for a preparation of the transmon in state 1 2) . b) Calculation 
based on Cavity-Bloch equations. 



wide range of detunings A rm when preparing the trans- 
mon in the |2) state and observing its decay into the |0) 
state. Three distinct maxima in the measured Q quadra- 
ture [Fig.[2ja)] located at the expected frequencies shifted 
by an amount s n from uo r are characteristic for the mea- 
surement of the n = 0, 1, 2 states of the transmon. The 
peaks appear successively in time, as the transmon se- 
quentially decays from |2) to |1) to the ground state |0). 
Sequential decay is expected due to the near harmonic- 
ity of the transmon qubit, for which only non- nearest- 
neighbor transitions are important [19]. The Q quadra- 
ture calculated from Cavity-Bloch equations is in good 
agreement with the measurement data and yields the en- 
ergy relaxation times of the first and second excited state 
Tl = 800 ± 50 ns and 7\ 2 = 700 ± 50 ns as the only fit 
parameters, see Fig.[2|b). The relaxation times are much 
longer than the typical time required to prepare the state 
using two consecutive 12 ns long pulses and allow for a 
maximum f-level population of 97%, limited by popula- 
tion decay during state preparation [32] . The relative 
difference between data and calculated transmission is 
at most 3% at any given point indicating our ability to 
populate and measure the f-level with high fidelity. 

To realize high-fidelity state tomography, arbitrary ro- 
tations in the Hilbert space with well defined phases and 
amplitudes are essential. Calibration of frequency, signal 
power and relative phases has to be performed based only 
on the population measurements of the qutrit states. To 
do so, we notice that the weak measurement partially 
projects the quantum state into one of its eigenstates 
|0), |1) or 1 2) in each preparation and measurement se- 
quence [2TJ[27 ? 28 . The average over many realizations of 
this sequence, which leads to the traces in Fig.[TJb), can 
therefore be described as a weighted sum over the contri- 
butions of the different measured states. This suggests 
the possibility of simultaneously extracting the popula- 
tions of all three levels from an averaged time-resolved 
measurement trace. Formally, the projective quantum 
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FIG. 3: (color online) Reconstructed transmon populations 
in a Ramsey oscillations experiment on the 1-2 transition us- 
ing a 5 MHz detuned drive field. For a given pulse delay |0) 
(blue dots), |1) (red diamonds) and |2) (green squares) popu- 
lations are extracted from a time resolved averaged measure- 
ment trace. The lines are calculated using Bloch equations. 



non-demolition measurement gives rise to the following 
operator, which is diagonal in the three-level basis and 
linear in the population of the different states at all times, 

Mj(t) = 5o(t)|0)(0|+5i(t)|l)(l|+5 2 (t)|2)(2|. (2) 

Here, a n (t) are the averaged transmitted field amplitudes 
for the states n sketched in Fig[TJa). The transmitted in- 
phase quadrature 

(I p (t)) = TV [pMj(t)] = poao (t)+Piai(£)+p 2 a 2 (£), (3) 

can be calculated for an arbitrary input state with den- 
sity matrix p and populations pi. Since any measured 
response is a linear combination of the known pure |0), 
|1) and 1 2) state responses weighted by pi, the popula- 
tions can be reconstructed using an ordinary least squares 
linear regression analysis, which pseudo-inverts Eq. Q 
for each time step U. The reconstructed populations 
show larger statistical fluctuations than in the two-level 
case [27] due to the pseudo-inversion of the ill-conditioned 
matrix used to calculate the pi from Eq. The statis- 
tical error is influenced by the distinguishability between 
the different traces, see Fig.[]Jb), and is minimized by op- 
timizing the measurement detuning. In contrast to full 
quantum state tomography (see below), this method does 
not require any additional pulses after the state prepa- 
ration. It is therefore used to find the U12 transition 
frequency and the pulse amplitudes needed to generate 
accurate pulses for tomography. 

The pulse amplitudes are extracted from Rabi- 
oscillations. To asses the precise value of CJ12, we perform 
a Ramsey experiment between the |1) and the |2) level, 
see Fig. [3] We apply a 7r-pulse at cjoi and then delay the 
time between two successive 7r/2 pulses applied at U12 be- 
fore starting the measurement. The theoretical lines are 
calculated based on a Bloch equation simulation with a 
dephasing time of the 2-level fitted to T| = 500 ns. 

Using quantum state tomography [TO] , the full den- 
sity matrix of the first three levels of a transmon can be 
reconstructed. This is achieved by performing a com- 
plete set of nine independent measurements after prepa- 
ration of a given state and calculating the density matrix 




FIG. 4: (color online) Measured Real and imaginary part of 
the reconstructed density matrices of |^ a ) = l/\/2(|l) — |2)) 
and \ySf b ) = l/>/3(|0> - |2)). The cyan cylinders indi- 

cate the standard deviations, typically 0.02. 



based on the measurements outcomes. Since the mea- 
surement basis is fixed by our hamiltonian Eq. 0, the 
state is rotated by applying the following pulses prior 
to measurement: I, (f )°\ (f )°\ (tt)°\ (§)* 2 , (f )* 2 , 

Ml 1 (f (f )^ Ml" to" where I denotes the 

identity and (6) l £ denotes a pulse of angle on the 
^/'-transition about the a-axis. For each of these uni- 
tary rotations (Uk) we measure the coefficients = 
Tr[pUkMiUl] by integrating the transmitted in-phase 
quadrature in Eq. Q over the measurement time [28] , 
i.e. implementing the measurement operator Mj = 
Jq Mi(t)dt. This relation is inverted to reconstruct 
the density matrix p by inserting the known operators 
UkMiU\. Note, that unlike in the preceding measure- 
ment of the populations only, we now extract a single 
quantity, (Ik), for each measured time trace. Quantum 
state tomography based on the simultaneous extraction 
of the populations of |0), |1) and |2) could potentially 
reduce the number of required measurements, but might 
come at the expense of larger statistical errors as dis- 
cussed above. 

Examples of measured density matrices are shown in 
Fig.Hfor the states \^ a ) = l/y/2(\l) - |2)) and \V b ) = 
1/Vo(|0) — |2)). A maximum likelihood estima- 

tion procedure has been implemented [33] . The extracted 
fidelities F = of 97 db 2% and 92 db 2%, respec- 

tively, demonstrate the high level of control and the good 
understanding of the readout of our three level system. 
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Considering the measured decay rates, the best achiev- 
able fidelity for the states \^ a /b) is 97 ± 1%. Preparing a 
set of 14 different states we measure an average fidelity 
of 96%, with a minimum of 92 ± 2% for the pure |2) 
state. The small remaining imperfections are likely due 
to phase errors in the DRAG pulses which affect both 
state preparation and tomography. 

We have demonstrated the preparation and tomo- 
graphic reconstruction of arbitrary three level states in 
a superconducting quantum circuit. Controlling and 
reading out higher excited states in these systems does 
broaden the prospects of using such circuits for future ex- 
periments in the domains of quantum information science 
and quantum optics. 

We thank Marek Pikulski for his contributions to the 
project. This work was supported by SNF, the EU 
projects SOLID and EuroSQIP and ETH Zurich. AB 
is supported by NSERC, the Alfred P. Sloan Foundation 
and is a CIFAR Scholar 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge Univertity Press, 
2000). 

[2] B. P. Lanyon, M. Barbieri, M. P. Almeida, T. Jennewein, 
T. C. Ralph, K. J. Resch, G. J. Pryde, J. L. O'Brien, 
A. Gilchrist, and A. G. White, Nature Physics 5, 134 
(2009). 

[3] M. Neeley, M. Ansmann, R. C. Bialczak, M. Hofheinz, 
E. Lucero, A. D. O'Connell, D. Sank, H. Wang, J. Wen- 
ner, A. N. Cleland, et al., Science 325, 722 (2009). 

[4] N. J. Cerf, M. Bourennane, A. Karlsson, and N. Gisin, 
Physical Review Letters 88, 127902 (2002). 

[5] T. Durt, D. Kaszlikowski, J. L. Chen, and L. C. Kwek, 
Physical Review A 69, 032313 (2004). 

[6] D. Kaszlikowski, P. Gnacinski, M. Zukowski, W. Mik- 
laszewski, and A. Zeilinger, Physical Review Letters 85, 
4418 (2000). 

[7] R. Inoue, T. Yonehara, Y. Miyamoto, M. Koashi, and 
M. Kozuma, Physical Review Letters 103, 110503 (2009). 
[8] A. Mair, A. Vaziri, G. Weihs, and A. Zeilinger, Nature 

412, 313 (2001). 
[9] G. Molina- Terriza, A. Vaziri, J. Rehacek, Z. Hradil, and 
A. Zeilinger, Physical Review Letters 92, 167903 (2004). 

[10] R. T. Thew, K. Nemoto, A. G. White, and W. J. Munro, 
Physical Review A 66, 012303 (2002). 

[11] G. Vallone, E. Pomarico, F. De Martini, P. Mataloni, and 
M. Barbieri, Physical Review A 76, 012319 (2007). 

[12] J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, 
Physical Review Letters 89, 117901 (2002). 

[13] K. B. Cooper, M. Steffen, R. McDermott, R. W. Sim- 
monds, S. Oh, D. A. Hite, D. P. Pappas, and J. M. Mar- 
tinis, Physical Review Letters 93, 180401 (2004). 

[14] E. Lucero, M. Hofheinz, M. Ansmann, R. C. Bialczak, 



N. Katz, M. Neeley, A. D. O'Connell, H. Wang, A. N. 

Cleland, and J. M. Martinis, Physical Review Letters 

100, 247001 (2008). 
[15] K. V. R. M. Murali, Z. Dutton, W. D. Oliver, D. S. 

Crankshaw, and T. P. Orlando, Physical Review Letters 

93, 087003 (2004). 
[16] J. Claudon, F. Balestro, F. W. J. Hekking, and O. Buis- 

son, Physical Review Letters 93, 187003 (2004). 
[17] S. K. Dutta, F. W. Strauch, R. M. Lewis, K. Mitra, 

H. Paik, T. A. Palomaki, E. Tiesinga, J. R. Anderson, 

A. J. Dragt, C. J. Lobb, et al., Physical Review B 78, 
104510 (2008). 

[18] A. Ferron and D. Dommguez, Physical Review B 81, 
104505 (2010). 

[19] J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. I. 
Schuster, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, 
and R. J. Schoelkopf, Physical Review A 76, 042319 
(2007). 

[20] J. A. Schreier, A. A. Houck, J. Koch, D. I. Schuster, 

B. R. Johnson, J. M. Chow, J. M. Gambetta, J. Majer, 
L. Frunzio, M. H. Devoret, et al., Physical Review B 77, 
180502(R) (2008). 

[21] A. Blais, R. S. Huang, A. Wallraff, S. M. Girvin, and 
R. J. Schoelkopf, Physical Review A 69, 062320 (2004). 

[22] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R. S. 
Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. 
Schoelkopf, Nature 431, 162 (2004). 

[23] M. Baur, S. Filipp, R. Bianchetti, J. M. Fink, M. Goppl, 
L. Steffen, P. J. Leek, A. Blais, and A. Wallraff, Physical 
Review Letters 102, 243602 (2009). 

[24] M. A. Sillanpaa, J. Li, K. Cicak, F. Altomare, J. I. Park, 
r. W. Simmonds, G. S. Paraoanu, and P. J. Hakonen, 
Physical Review Letters 103, 193601 (2009). 

[25] L. DiCarlo, J. M. Chow, J. M. Gambetta, L. S. Bishop, 
B. R. Johnson, D. I. Schuster, J. Majer, A. Blais, L. Frun- 
zio, S. M. Girvin, et al., Nature 460, 240 (2009). 

[26] A. A. Abdumalikov Jr., O. Astafiev, A. M. Zagoskin, 
Y. A. Pashkin, Y. Nakamura, and J. S. Tsai, 
arXiv: 1004.2306 (2010). 

[27] R. Bianchetti, S. Filipp, M. Baur, J. M. Fink, M. Goppl, 
P. J. Leek, L. Steffen, A. Blais, and A. Wallraff, Physical 
Review A 80, 043840 (2009). 

[28] S. Filipp, P. Maurer, P. J. Leek, M. Baur, R. Bianchetti, 
J. M. Fink, M. Goppl, L. Steffen, J. M. Gambetta, 
A. Blais, et al., Physical Review Letters 102, 200402 
(2009). 

[29] P. J. Leek, M. Baur, J. M. Fink, R. Bianchetti, L. Steffen, 
S. Filipp, and A. Wallraff, Physical Review Letters 104, 
100504 (2010). 

[30] F. Motzoi, J. M. Gambetta, P. Rebentrost, and F. K. 

Wilhelm, Physical Review Letters 103, 110501 (2009). 
[31] F. Motzoi, J. Gambetta, and F. Wilhelm, in preparation. 
[32] J. M. Chow, J. M. Gambetta, L. Tornberg, J. Koch, L. S. 

Bishop, A. A. Houck, B. R. Johnson, L. Frunzio, S. M. 

Girvin, and R. J. Schoelkopf, Physical Review Letters 

102, 090502 (2009). 
[33] D. F. V. James, P. G. Kwiat, W. J. Munro, and A. G. 

White, Physical Review A 64, 052312 (2001). 



